REMARKS ON A QUASI-LINEAR MODEL OF THE 
NAVIER-STOKES EQUATIONS 
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Abstract. Dinaburg and Sinai recently proposed a quasi-linear model of 
the Navier-Stokes equations. Their model assumes that nonlocal interactions 
in Fourier space are dominant, contrary to the Kolmogorov turbulence phe- 
nomenology where local interactions prevail. Their equation corresponds to 
the linear evolution of small scales on a background field with uniform gradi- 
ent, but the latter is defined as the linear superposition of all the small scale 
gradients at the origin. This is not self-consistent. 

Dinaburg and Sinai [J recently proposed a quasi-linear approximation of the 
Navier-Stokes equations which they feel preserves the basic character of the Navier- 
Stokes nonlinearity. They prove existence and uniqueness of solutions to their 
model, for special cases [2- Here, we show that their equation is identical to that 
governing the linear evolution of small scales on an infinitely large scale flow with 
uniform gradient (eqn. below) . This is a direct consequence of their assumption 
that nonlocal interactions in Fourier space dominate. That assumption is contrary 
to a large body of phenomenological, experimental and numerical studies of the 
turbulent energy cascade where local interactions are thought to dominate (see 
e.g. j3j). On that basis alone, it seems unlikely that the model would preserve 
the basic character of the Navier-Stokes nonlinearity, but the model also contains 
a basic inconsistency. It defines the velocity gradient of the large scale field as 
the net velocity gradient induced by all the small scales at the origin (eqn. 1)12(1 
below). This is not self-consistent since small scales do not uniformly distort larger 
scales. A self-consistent mean field theory, where the large scale flow results from 
the nonlinear interactions of the small scale fluctuations, instead of their linear 
superposition, is not possible in the nonlocal limit. 

Consider the Navier-Stokes equations for incompressible flow of a viscous fluid 
in three-dimensional space M'^: 

du , ^. ^ 
— + {u-S/)u + Vp^vV-u, 
(1) dt 

V • u = 0, 

where u = u(x, t) — (mi(x, i), U2(x, t),u^(x, t)) is the velocity vector field, p — p{x, t) 
is the kinematic pressure, x — (xi,X2,x^) is the position vector and > is the 
kinematic viscosity. 

Let u = A-x + v where A ~ A[t) is a matrix in x such that dA/dt + A- A 
is symmetric and tr(A) = in order for A • cc to be a solution of the incompressible 
Navier-Stokes equations The base field A ■ x has uniform gradient and v = 
v{x,t) = {vi{x,t),V2{x,t),V3{x,t)) is a velocity perturbation. 
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Substituting u = A-x + vi'a^ and linearizing in v {i.e. omitting (w- V)u) yields 
dv 



(2) 



+ (^{A ■ x) ■ y^v + A ■ V + ^ vV 



or, in indicial notation, 



d^n . dvi ^ ^ ^ dp d^vi 



L/C (JJLjYi (JJLl C/ 777 L/ X 7] 



0. 



dxi 

We will use the convention of summation over repeated indices l,m,n = (1,2,3). 
The self-advection of the background field is the gradient of a potential that has been 
absorbed into the pressure, since dA/dt + ^ • A is symmetric. Incompressibility of 
the base field requires tr(v4) = An = 0. An equation for the pressure can be derived, 
in a standard manner, by taking the divergence of the momentum equation in ^ 
and using incompressibility to obtain 

(4) V'P=^ = -2A..,*' 



OXldXl OXm 

Linear problems of the form ^ have been considered by many authors since 
Kelvin. There has been a renewed interest in such analyses in more recent times 
linked to a simple model of the elliptical instability (see ^4, ^-nd references therein). 

Kelvin modes. Kelvin noticed that a generalized Fourier analysis could be 
used to solve equation |5J). He proposed to look for solutions of the form 

(5) w(a;,i) ={;(i)e''=(*)-" 

i.e. Fourier modes with time-dependent wave- vectors fc £ M"^. The motivation 
for this ansatz is that a Fourier mode initial condition proportional to e^^'^ is 
rotated and stretched uniformly by the background field with uniform gradients. 
Hence, it remains in the form of a Fourier mode, albeit with an evolving wavevector 
k. Substituting © together with p{x,t) = p{t)e'^^'^^''^ into ©, and using (gj) to 
eliminate the pressure, leads to the coupled ordinary differential equations: 

dk , , 

(6) ^ ^ = 

(7) ^^=.^\k\^i,^A-v + 2^{k-A-v). 

The incompressibility constraint V • w = requires k{t) ■ v{t) — Q. This is satisfied 
automatically since the pressure is determined from igj, provided that the initial 
conditions are such that fc(0) • w(0) = 0. The most interesting solutions of these 
equations, perhaps, occur when the base flow A ■ x has closed streamlines. Then 
k{t) is oscillatory and v{t) can grow exponentially through a parametric instability 
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Fourier transform. Alternatively, one can proceed with a direct Fourier trans- 
form of the equations. Let v{k,t) be the Fourier transform of the vector 

(8) v{x,t) = I v{k,t)e'^-''dk, 
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then 



(9) A^n^n 1^ — / AjnnTr, — ( kmVlikjt)^ e'''' ^ 

dx,n Jrs Okn V / 

which, in our case, simpUfies to 



dk, 



OXm JR3 OKn 

because Amm = 0. It foUows that the Fourier transform of equation ||2J) reads 
, , dvi{k,t) dvi 17 |2- /I - o ^' /; /I - \ 

(11) 'q^ kjyiAjnn = —v\k\ Vi — AijjiVm + ( "i ) ' 

where the pressure and the incompressibihty constraint have been eliminated using 
provided the initial conditions satisfy kivi{k,0) = 0, Vfc. Solving equation l|ll|) 
by the method of characteristics, we directly recover Kelvin modes and equations 

®, 0- 

Equation (|ll|l is identical to Dinaburg and Sinai's equation (11) except for a 

different definition of the matrix A, theirs is minus the transpose of ours: A'^'^' — 
—A^. Note that equation 1)11(1 has been obtained without imposing the point 
symmetry v{x,t) = —v{—x,t), assumed by Dinaburg and Sinai, although the base 
flow A ■ X does satisfy that symmetry. The point symmetry implies that v{k, t) = 
—v{—k,t). Now, v{k,t) — v{—k,t) since v{x,t) is real, where the overline denotes 
complex conjugate. Hence, v{k,t) must be pure imaginary if the point symmetry 
is imposed. So Dinaburg and Sinai's v{k, t) is related to our w(fc, t) as v^^^^ (fc, t) = 
—iv{k,t). 

Dinaburg and Sinai define 

(12) A^Z'^ =- f ^kMk, t)dk = - 

J OXm 

while, if we write our base flow as U{x, t) — A{t) ■ x, then — A^nXn and 

OXji 

This explains the differences in the definitions of the velocity gradient A. It also 
points to a basic inconsistency of the Dinaburg-Sinai model. Equation (|ll|l corre- 
sponds to the distortion of small scales by an infinitely large scale flow with uniform 
gradient, but Dinaburg and Sinai define the gradient of the infinitely large scale 
velocity field as the local gradient at a; = resulting from the linear superposi- 
tion of all the small scale gradients. Considering two distinct Fourier modes with 
wavevectors A:*^^^ and k'^^\ for instance, with Ifc*^^^! < \k^^^\, it does not make sense 
to have the small scale fc^^^ participating in the uniform large scale distortion of 
the larger scale k'--^K Furthermore, the mode /c^^-* does not self-distort because of 
the incompressibihty constraint (so (v ■ V)u = for a single Fourier mode). Hence, 
only larger scales should contribute to the approximation of uniform distortion 
of a given small scale. In other words, only \k'\ <^ \k\ should contribute to the 
large scale gradient distorting the mode with Fourier wavevector k. This can be 
seen also at a technical level in the derivation of the Dinaburg-Sinai model. In 
their treatment of the convolution integral representing the Fourier transform of 
the Navier-Stokes nonlinearity, Dinaburg and Sinai make the assumption that the 
integral is dominated by highly nonlocal interactions, i.e. by the domains \k'\ <C |fc| 



4 



FABIAN WALEFFE 



and — fc'l ^ \k\. Considering the domain \k'\ ^ |A:|, they make the following type 
of approximation, for instance, 



(see eqn. (6) in J^). The domain of integration for the integral on the right-hand 
side should be restricted to the ball ~ {k' G M"^ : |fc'| < e}, with e <C |fc|, for 
self-consistency, but Dinaburg and Sinai do not specify the domain of integration. 
In their later studies of finite dimensional approximations, they sum over all modes, 
in other words, they integrate over M'^ instead of B^. If the integral, and the finite 
dimensional sums, were correctly restricted, the model would not be closed, or the 
infinitely large scale flow would have to be specified and the model would become 
linear and identical to Hll() . with A{t) specified independently from the small scales. 

It is natural to wonder whether one could replace the Dinaburg-Sinai model by a 
mean field theory where the large scale fiow results from the nonlinear interactions 
of the small scales. The base field U — A ■ x is a. very singular fc = mode whose 
generalized Fourier transform U{k) — i{A ■ V/c)(5(fc). Here 6{k) = S{ki)6{k2)S{k3) 
is a product of Dirac delta functions and is the gradient operator in fc-space 
with d6{k)/dki — S' {ki)6{k2)6{k3), where S'{ki) is the generalized derivative of 
the delta function, and similarly for derivatives with respect to ^2 and k^. The 
only nonlinear interactions that can create a fc = mode consist of any mode fc' 
interacting with its complex conjugate — fc', but such interactions vanish because 
of incompressibility. For instance, the convolution integral on the left-hand side of 
(PHl when fc = 0, 



for any regular u(fc') because fc^u„(— fc') = from incompressibility. The vanishing 
of such interactions is also related to Galilean invariance since a fc = mode could 
also correspond to a constant velocity. 
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